Single-copy entanglement in critical quantum spin chains 
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We introduce the single-copy entanglement as a quantity to assess quantum correlations in the ground state 
in quantum many-body systems. We show for a large class of models that already on the level of single speci- 
mens of spin chains, criticality is accompanied with the possibility of distilling a maximally entangled state of 
arbitrary dimension from a sufficiently large block deterministically, with local operations and classical commu- 
nication. These analytical results - which refine previous results on the divergence of block entropy as the rate 
at which EPR pairs can be distilled from many identically prepared chains, and which apply to single systems 
as encountered in actual experimental situations - are made quantitative for general isotropic translationally 
invariant spin chains that can be mapped onto a quasi-free fermionic system, and for the anisotropic XY model. 
For the XX model, we provide the asymptotic scaling of ~ (1/6) log 2 (L), and contrast it with the block en- 
tropy. The role of superselection rules on single-copy entanglement in systems consisting of indistinguishable 
particles is emphasized. 

PACS numbers: 03.75.Ss, 03.75.Lm, 03.75.Kk 



Quantum phase transitions of second order are accompa- 
nied with a divergent length scale: this is the classical cor- 
relation length, the characteristic length associated with the 
two-point correlation function 1,1]. Recently, it has increas- 
ingly become clear that one should expect additional insight 
in the scaling of quantum correlations present in the ground 
state of a many-body system at or close to a quantum phase 
transition by expressing them in terms of entanglement prop- 
erties HElSS|lQ,S0[l3[li. Entan g lement ' af- 
ter all, plays a fundamental role in quantum phase transi- 
tions at zero temperature. The theory of entanglement in turn 
- developed in the quantum information context - provides 
tools to characterize and quantify genuine quantum correla- 
tions in contrast to correlations that occur in states that can 
be prepared with mere local preparations together with clas- 
sical communication (LOCC). In particular, one finds that in 
one-dimensional non-critical harmonic 0. 0. 01 or quantum 
spin systems 0, 01, the degree of entanglement of 
a block of L systems, quantified in terms of the entropy of 
the reduction, typically saturates for large block size, with 
higher-dimensional "entropy-area laws" |4J|. In contrast, in 
critical spin systems or in fermionic systems, the entropy of a 
reduction has logarithmic corrections as L — > oo 
These findings are consistent with expectations from confor- 
mal field theory |2]. Such a behavior of the block entropy has 
also been related to the performance of DMRG simulations 
of ground state properties. This von-Neumann entropy of a 
block quantifies the rate at which one can asymptotically dis- 
till maximally entangled qubit pairs under LOCC, when one 
has infinitely many identically prepared many-body systems 
at hand full. 

Yet, in several contexts, in particular for condensed-matter 
systems, this asymptotic notion of entanglement implicitly re- 
ferring to joint operations on many identical systems may not 
always be the most appropriate one. Instead, one may ask: 
does a single specimen of a critical infinite system already 
contain an infinite amount of entanglement? This will be the 



central question of this paper. We introduce the single-copy 
entanglement to quantify the quantum correlations in criti- 
cal and non-critical many-body systems. More specifically, 
compared to the divergence of the block entropy, we ask the 
stronger question whether a single spin chain already contains 
an arbitrary amount of entanglement, such that from a single 
specimen a maximally entangled state of arbitrary dimension 
can be distilled. 

We will make the argument quantitative by analytically 
considering a general framework of translationally invariant 
quantum spin models. As examples in which criticality is in 
one-to-one correspondence with a divergent single-copy en- 
tanglement, we consider isotropic spin models, as well as the 
XY-model. For the isotropic XY model we establish the exact 
asymptotic scaling behavior of ~ (1/6) log 2 (L), and relate it 
to the block entropy 00]. The results can also be conceived 
as statements concerning the divergence of fine-grained en- 
tanglement fl2ll . 

Single-copy entanglement. - Let us consider a one- 
dimensional quantum spin system, associated with a Hilbert 
space Ti — ((D 2 )®", with a translationally invariant Hamilto- 
nian. We distinguish a block of length L of consecutive sys- 
tems of the chain. So we have a bi-partioning n\L, the whole 
system being in a pure state p = 

There are several meaninful definitions of single-copy en- 
tanglement. We will primarily be concerned with the ques- 
tion: running a physical device once, a maximally entangled 
state of what dimension can be distilled from a single speci- 
men with certainty? Hence, the state has a single-copy entan- 
glement Ei (p) = log 2 (M), with respect to the bi-partitioning 
n\L, if p can be deterministically transformed under LOCC 
into IV'mXV'mI' i- e -' a maximally entangled state with state 
vector |V>m) = (|1, 1) + ... + \M, M))/VM, so if 

p — ► |^>M)(V'Af| under LOCC. (1) 

This is the non-asymptotic analogue of the entropy of entan- 
glement of the reduction associated with a block of length L. 
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Denote with a[, a^ L the non-increasingly ordered eigen- 
values of the reduced state with respect to a block of length 

L, then Eq. © holds true if and only if Q J2k=i a k — 
K/M for all 1 < K < M, so obviously, if and only if 
a\ < 1 /M. In other words, the transformation is possible if 
the reduction is more mixed in the sense of majorization than 
the reduction of the maximally entangled state of dimension 
M x M. Given a[, the single-copy entanglement is nothing 
but Ei(p) = \og 2 ([(a{)~ 1 \ ). A variant is one allowing for 
probabilistic protocols. For a state p we say that 



E p (p) = sup^p fe log 2 (A/ fc ), 

k=0 

such that p can be transformed under LOCC into the ensem- 
ble {(p k , |-0M fc )(-0A/J) : k = 0, 1, ...}. This is the average 
entanglement that can be distilled, allowing for maximally en- 
tangled states of different dimension with certain probabilities 
111411 . This rate is then the solution of a linear program 
By definition, we have that E\ (p) < E p (p) < S(ti n \ L [p\), 
where the last inequality follows from the fact that the entropy 
of a reduction bounds the rate of any (asymptotic) distillation 
protocol. 

Finally, note that for the single-copy entanglement, super- 
selection rules (SSR) play a crucial role, notably a SSR with 
respect to particle number conservation. In the presence of 
SSR, Ef SR has to be understood as referring to a probabilistic 
transformation distilling entanglement-EPR states with state 
vector |^i SSR ) = (|0,1)|1,0) + |l,0)|0,l))/ v / 2underLOCC 
and SSR. We will now consider the behavior of the single- 
copy entanglement in the limit of large L for critical and non- 
critical spin chains. 

Single-copy entanglement in general quantum spin chains. 
- We start from the general set of translationally invariant 
quantum spin systems that is as in Ref. 0] mapped onto 
a fermionic quadratic Hamiltonian under a Jordan- Wigner 
transformation. This model embodies a large class of spin 
models, including the anisotropic and isotropic XY-models as 
important special cases. Hence, the Hamiltonian is 



ff =EE 

g=0,l kj=l 



Bj-k Aj- 



n 

i<k + q 
l<j+l-q 
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where eh, vf. denote the Pauli operators associated with 
site k = 1, n, equivalent with the fermionic Hamiltonian 



n 

= ^2 [o]Aj_fcOfc + a\Bj_ k a\ - ajBj- k ak 



The fermionic operators obey {aj,cik} — and {a],dfc} = 
Sj t k- The Hamiltonians are related via a Jordan- Wigner 
transformation leading to the Hermitian Majorana operators, 

msi-i = (Uj<i ^ ™2i = {Uj<i where 

a,j = (rhzj-i — im,2j)/2. Translational invariance, peri- 
odic boundary conditions, and Hermicity are inherited by 
Aj,Bj £ 1R satisfying Aj = A_j, and Bj = —B_j for 



j = 1 — n, n — 1. For simplicity, we assume that there 
exists aw e N such that Aj = B 3 ■ = for j > w. This 
model will be our starting point. For all isotropic instances, 
and also for the full XY model we will be able to identify 
when the single-copy entanglement is indeed logarithmically 
divergent. Then, one may distill a maximally entangled state 
of any dimension from a single specimen of the chain with 
certainty, containing in this sense an "infinite single-copy en- 
tanglement" 12511 - We will make use of the powerful methods 
of Toeplitz determinants Q |^ fT3l . This path is yet in our in- 
stance complicated by the fact that we do not only consider 
isotropic models, and that in contrast to the block entropy the 
largest eigenvalue cannot straightforwardly be expressed as an 
integral of a Toeplitz determinant. The starting point, yet, is 
the familiar one for assessing spin systems: The ground state 
of this system is a fermionic Gaussian, i.e., quasi-free, state 
and is completely specified by the second moments of the 
Majorana operators. These operators satisfy rhj = rhj and 
{rhj, rhk} = 2(5^. The second moments can be collected in 
a correlation matrix 7 6 j^2nx2n^ trfpmjm^] = S j & + ijj.k- 
This matrix is skew-symmetric. The entanglement properties 
of the block of length L can now be inferred from a principal 
submatrix 7^ e j^2Lx2L q £ ^ c orrel a.tion matrix 7. We con- 
sider the entries of jl in the limit of an infinite chain n — > 00. 
Then, 72, is a block Toeplitz matrix, the Z-th row, I = 1, L, 
being given by (M ; _i, M z _ 2 , M , M Z _ L ), with 2x2- 
blocks Ml -1 , ■ ■ ■ , Mi_£ that are found to be 
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For no anisotropy, i.e., Bj = for all j, this matrix is a tensor 
product of a symmetric matrix and a unit skew-symmetric one. 
In generality, we have for this model, g(k) := A(fc)/|A(fc)|, 
A(fc) := A +2 £J =1 Aj cos(jk)-4i EJ =1 Bj sin(jfc). This 
matrix jl can be brought into a standard normal form of 
a skew-symmetric matrix with an O E 0(2L) preserving the 
anticommutation relations, 



e 

1=1 



pi 
-Pi 



This defines the quantities pi, p^ £ [0,1]. Such normal 
mode decompositions have been employed both to evaluate 
correlation functions fl8ll and the block entropy 0, 0]- 
From now on we will be concerned with the largest eigen- 
value of the reduction of a block of length L. All eigenval- 
ues a[, a| L of the reduction are given by {a{, a^} = 
{Y[h=i (1 i / i /)/2}- We will be looking at the behavior of the 
largest eigenvalue a\ for large L. This largest eigenvalue is 
given by aj = U^Lii 1 / 2 + Vi/ 2 )' or 

a{ = det[(l£ + |T £ |)/2], 

\ t l\ = (TZTl) 1 / 2 , where T L is the L x L Toeplitz matrix, 

with l-th row being given by ? ^ ^+2? ■■■;^0; ■•*) ^L— i)* 

The numbers pi, ...,Pl are the singular values of T £ . This 
matrix T £ , satisfying |T £ | < 1^, is generally not symmetric, 
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as a consequence of the anisotropy of the model. Moreover, in 
contrast to the matrix Tl itself, Ij, + \Tl\ is not Toeplitz. In 
order to show that the single-copy entanglement is logarithmi- 
cally divergent, we will make use of appropriate bounds that 
retain this property: whenever the Aq, A w , Bq, B w are 
such that one can prove that the sequence of L x i-Toeplitz 
matrices Tl satisfies 



■iog|det[r L ]| = n(io g (i)) 



(2) 



(using Landau notation 12 CUP using a Fisher-Hartwig- 
statement men one can indeed conclude that E\ = 

fi(log(L)), i.e, the single-copy entanglement diverges at least 
logarithmically with increasing block length L. This follows 
from the following chain, 

- logdetpx + \T L \)/2] > -ilogdet[(l L + TjT L )/2] 
> - 1 log dot [TlT L ] = -\ log | dot [T L ] \ 



12711 . where we also have made use of the concavity of the 
logarithm. So, whenever Eq. (0 holds, for appropriate length 
of the block L, a maximally entangled pair of any dimension 
can be distilled from a single specimen of the spin chain. 

Isotropic models. - This case of Bq, B w = is partic- 
ularly transparent. Here, the asymptotics in L of the deter- 
minants det^fj,^] of the L x L-Toeplitz matrices M x l '■= 
ixt + (1 — x 2 ) 1 / 2 Tl is known for all x e (0, 1), using a 
Fisher-Hartwig statement. This small detour to infer about 
— log | det [Tl] \ - corresponding to the case x = - is needed 
as the Fisher-Hartwig-conjecture has not been proven yet for 
this case. In general, one can identify the asymptotic behavior 
of determinants of Toeplitz matrices by investigating the so- 
called symbol, see footnote Bill . The symbol associated with 
the Toeplitz matrices M x< l is given by 

G x {k)=ix + {l-x 2 ) 1 ' 2 g(k), 

with g as defined above. For this class of isotropic models, an 
explicit factorization of the symbol is known H, see footnote 
12211 . It follows hence from proven instances of the Fisher- 
Hartwig conjecture that there exists a c > and an xq <G (0, 1) 
such that ^ 



log | det[M x , £ ]| = c x log(L) + o(log(L)) 

with c x > c for all x <G (0, xq), whenever the function g 
is discontinuous in [0, 2ir], where the jumps reflect the Fermi 
surface. From this - and using that Tl has real eigenvalues - it 
follows that the system has a logarithmically divergent single- 
copy entanglement if the system is critical 12311 . For example, 
for the XX model this analysis immediately delivers a loga- 
rithmically divergent single-copy entanglement, whenever the 
system is critical. 

Anisotropic XY-model. - For the XY-model we can con- 
clude that the single-copy entanglement is logarithmically di- 
vergent if and only if the system is critical. For this model, 
we have that Aq = — l,Ai = a/2 and B-\ = — B\ 



~.n 1, and elsewhere. For 7 = 0, we obtain the XX- 
model (the isotropic XY-model), for a = 1, 7 = 1 the crit- 
ical Ising model. Along the line 7 £ [—1, 1], a = 1 the 
anisotropic model is critical. Then, we encounter a generally 
non-symmetric matrix Tl- The associated symbol is given by 



a cos(fc) — 1 + iaj sin(fc) 
((acos(fc) - l) 2 + 7 2 a 2 sin 2 ^)) 1 ^ 



For 7^0 and 1/a G (0, 1), the symbol is continuous, and one 
finds a saturating block entropy Q] (and hence a saturating 
single-copy entanglement). Along the critical line a = 1, 7 € 
(—1, 1), in turn, we can identify the explicit factorization of 
the discontinuous symbol. There is a single discontinuity at 
ki = I22I, and in the terms of footnote il we find /3i = 
1/2, so that g(k) can be decomposed as 

g(k) = <f>(k)t 1/2 (k), 

where is a continuously differentiable function. For the 
case of a single discontinuity and ct\ =0, the Fisher- 
Hartwig conjecture has been proven for any /3i g C with 
5f ) < 5/2 lll9ll , including our case at hand. Hence, we find 
-log|det[T L ]| = fi(log(L)), and hence E x = n(\og(L)). 
Together with the result of the subsequent section this shows 
that the single-copy entanglement of the XY-model is loga- 
rithmically divergent exactly if the model is critical. Note that 
this implies also a less technical alternative proof of the log- 
arithmic divergence of the block entropy in the critical XY 
model. 

Scaling of single-copy entanglement in the XX-model. - In 
the light of these findings, it is interesting to see how the exact 
asymptotic behavior is compared to that of the block entropy, 
including prefactors. We make this specific for the isotropic 
XX-model, where now Tl = Tj. The technicality when eval- 
uating a} = det[(l£ + |Tl|)/2] that we encounter here is 
that the function /:©-►€, f(x) := log 2 ((l + |ar|)/2), 
is not analytic. So before we can exploit Fisher-Hartwig-type 
results, we have to approximate a{ with sequences based on 
functions with appropriate continuity properties. We can take 
any functions /» : (D x R + — > (D which are analytic on 
{z 6 <D : 3(2) < 6} for a S > 0, such that on the real 
axis lim^o f*(x, S) = f(x, 0) for xeE. Take, e.g., 

f„{z,S) :=log(l/2+(z 2 + ( 5 2 ) 1 / 2 /2). 

We are then in the position to identify the asymptotic behav- 
ior of the single-copy entanglement. This can be done simi- 
larly to Ref. [7] using the characteristic polynomial F : (D — > 
€ of T L defined as F(\) := dct[Al L - T L \. the function F 
is meromorphic, and all zeros are in the interval [—1,1]. One 
can hence write 



lim lim 

J\o £ \o 2ni 



dzMz,6) 



F'{z) 
F(z) 



(3) 



where the integration path is chosen to enclose the interval 
[—1, 1], with path from (—1 — S + is, 1 + S + ie), towards 
the negative real numbers along a circle segment with radius 
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(5/2, then (1 + S — ie, — 1 — 5 — ie), and again along a circle 
segment to —1 — 5 + ie, such that lira^o = d. The symbol 
of Al — Tl with factorization as in Eq. @ for the XX-model 
is known @], see footnote l22ll . Using a Fisher-Hartwig state- 
ment, we find that the linear terms in L do not contribute, us- 
ing Cauchy's theorem and using that lim^o /»(±1, S) = 0, 
and finally arrive at 

d = log(i 4 /' ^[a + N)/^ +o( , OB(L)) . 

7r j_ 1 1 a; 

This in turn finally implies that whenever 1/a G [—1,1] and 
the XX-model is critical, we observe the scaling behavior 

Ei = ilog 2 (i)+o(log(L)), 

independent of a; it saturates in the non-critical case. This 
result is astonishing: the single-copy entanglement does not 
only diverge, but has up to a factor of two the same asymp- 
totic behavior as the entropy of entanglement scaling as S = 
(1/3) log 2 (L) + o(log(X)). Half of the asymptotically distil- 
lable entanglement is hence already available on the single- 
shot level. 

Outlook and summary. - Finally, let us comment on the 
crucial role of SSR for the single-copy entanglement. This is 



relevant, e.g., when assessing the single-copy entanglement in 
the hard-core limit of the Bose-Hubbard model (infinite repul- 
sion energy) There, the Hamiltonian is isomorphic to the 
XX-model, via the mapping <rj = bj + bj, &^ = —i(bj — St), 

and &j = 1 — 2bjbj for each site j. Yet, the concept of entan- 
glement is different due to the presence of a particle number 
conservation SSR in the former case: Transformations under 
LOCC have to be replaced by those under LOCC+SSR. The 
single-copy entanglement in the above sense can however still 
be efficiently evaluated in L; and these superselection rules 
must be respected when assessing single-copy entanglement. 

In this paper, we have fleshed out the notion of single-copy 
entanglement in quantum spin chains. Such a notion is the 
appropriate one when one is not interested in the entangle- 
ment properties of an asymptotic supply of a identically pre- 
pared many-body systems, but of single specimens. It is the 
hope that these findings also serve as a guideline when assess- 
ing entanglement in actual experimental situations, let it be in 
condensed-matter systems or in systems of ultracold atoms in 
optical lattices. 
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